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We propose a microscopic approach describing the interaction of an ideal gas of hydrogenlike 
atoms with a weak electromagnetic field. This approach is based on the Green-function formalism 
and an approximate formulation of the method of second quantization for quantum many-particle 
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and damping rate of electromagnetic pulses on the microscopic characteristics of the system are 
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I. INTRODUCTION 

The possibility of observation of ultraslow light pulses 
in vapors of alkali-metal atoms at extremely low temper- 
atures was shown in 1999 [l|. To the present moment, 
there are quite a lot of theoretical works devoted to the 
description of this phenomenon (probably the earliest of 
which is represented in Ref. |2|]). Naturally, for the most 
complete description of electromagnetic waves slowing in 
a Bose-Einstein condensate (BEC), the development of 
an appropriate microscopic approach is needed. 

But, while constructing this approach, some problems 
appear. First of all, one needs to use the method of 
second quantization for a microscopic description of the 
response of the BEC to the perturbation by the external 
electromagnetic field. The second-quantization method 
uses the creation and annihilation operators of atoms. 
This fact requires consideration of atoms as elementary 
objects, which do not have internal structure. However, 
atoms are compound objects with a rather complicated 
internal structure. Moreover, in the processes of interac- 
tion of neutral atoms with an electromagnetic field, this 
structure plays a crucial role. In particular, it leads also 
to a sufficient reduction of the group velocity of propa- 
gation of electromagnetic waves in a BEC. 

Thus, a problem appears concerned with the introduc- 
tion of the creation and annihilation operators, which 
must preserve the information about the internal struc- 
ture and, first of all, about the energy state of the atoms. 
In other words, we discuss the formulation of the second- 
quantization method for a many-particle system in the 
presence of bound states of particles (atoms). This prob- 
lem is solved in Ref. [||. In that work, an approximate 
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formulation of the second quantization method for such 
a system was constructed in the case when the binding 
energy of compound particles is much larger than the 
kinetic energy. The possibility of development of such 
a theory is demonstrated by considering a system that 
consists of fermions of two different kinds and bound 
states (atoms) that can be formed only by two fermions 
of different kinds (ions and electrons). The choice of this 
model in Ref. Q was dictated not by methodological dif- 
ficulties, but to get more visual results. In that work, 
the Hamiltonians of interaction for particles of different 
kinds (including bound states) were constructed. The 
Hamiltonian of interaction of all particles in the system 
(including neutral compound particles) with the electro- 
magnetic field was found. Because of this, the authors 
constructed a nonrelativistic quantum electrodynamics 
of systems which may consist of bound states of parti- 
cles. Within this theory, the Maxwell-Lorentz equations 
were obtained. These allow study of the response of a 
system consisting of charged "elementary" particles and 
bound states to perturbation by an external electromag- 
netic field. For solving this problem it is natural to use 
the Green-function formalism (see, e.g., Refs. 0, [f|). 

In Ref. [j| the Green-function method was generalized 
to the case of description of this response of the sys- 
tems mentioned above. In terms of Green functions, the 
expressions for such macroscopic characteristics as con- 
ductivity, permittivity, and magnetic permeability of a 
hydrogenlike plasma were found. 



In the present paper, the formulas that were found in 
Refs. 0, |j| are used to study the phenomenon of elec- 
tromagnetic pulses slowing in a BEC. According to the 
main results from Refs. p, 0], let us briefly summarize 
some basic principles of the proposed approach. 
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II. HAMILTONIAN OF HYDROGENLIKE 
ATOMIC GAS IN AN EXTERNAL 
ELECTROMAGNETIC FIELD 

Usually, BEC is studied in gases of alkali- metal atoms. 
It is known that the internal structure of alkali metals is 
similar to the internal structure of a hydrogen atom (hy- 
drogenlike atoms). Thus, such atoms may be considered 
as bound states of two elementary particles (the valence 
electron and the atomic core). In this case, the results 
of Refs. [H, Q become suitable for a description of the 
interaction of a weak electromagnetic field with a BEC. 
Let us note that in Ref. [f| the Green functions of hydro- 
genlike plasma were found. In other words, the contri- 
bution of free fermions and bound states (atoms) to the 
Green functions was considered. But the recent theoret- 
ical investigations Q of a low-temperature hydrogenlike 
plasma in the equilibrium state show that the density of 
free fermions in such a system is exponentially small in 
comparison with the atomic density. Therefore, at low 
temperatures one may neglect the free-fermion contribu- 
tion in the expressions for the Green functions. Taking 
into account this fact, the system Hamiltonian may be 
written as follows: 

n(t) = n + n int + v(t), H = n ph + H p , (i) 

where H. p h is the Hamiltonian for free photons and 7ii n t is 
the Hamiltonian of interaction between atoms. Note that 
we neglect the Hamiltonian H p h below. We also consider 
the system of hydrogenlike atoms as an ideal gas. Thus, 
we do not take into account the presence of the interac- 
tion Hamiltonian 7ii n t in the next calculations. Its mi- 
croscopic expression is found in Ref. Q in the framework 
of the developed formulation of the second-quantization 
method. This Hamiltonian is constructed on the basis of 
the Coulomb interaction between particles that can form 
a bound state (hydrogenlike atom). The explicit form 
for it in a low-energy approximation is given in the Ap- 
pendix. We use it to find the limits of applicability for 
the model of the ideal gas. 

The operator 7i p in Eq. JT]) is the Hamiltonian for free 
particles (atoms), 
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where iy a (x) and fj a (x) are the creation and annihilation 
operators of hydrogcnlikc (alkali-metal) atoms with the 
set of quantum numbers a at the point x, respectively; 
e a is the atomic energy corresponding to this state; and 
to is the atomic mass. 

The operator V(t) in Eq. fl} represents the Hamil- 
tonian that describes the interaction of atoms with an 
electromagnetic field, 

V(t) = -~ j dx A(x,i)j(x) + I dx V (x,t)ff(x), (3) 



where A(x, t) and </?(x, t) are the vector and scalar poten- 
tials of the external electromagnetic field, and j(x) and 
<t(x) are the operators of the current and charge densities 
for neutral atoms, respectively (see Ref. Q for details), 

*w = ^EE eiX(p '^ p) ^(p-p')^(p)^(p'), (4) 
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where V is the volume of the system. The matrix el- 
ements of the charge and current densities in Eqs. ((4]) 
and ([5]) can be expressed in terms of the atomic wave 
functions </? a (x): 
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where e is the electron charge absolute value, and m p 
and m e are the masses of the atomic core and electron, 
respectively (to = m p + m c ). 



III. RESPONSE OF A GAS OF 
HYDROGENLIKE ATOMS TO PERTURBATION 
BY AN EXTERNAL ELECTROMAGNETIC 
FIELD: GREEN FUNCTIONS 

Here we find the expressions that describe the influ- 
ence of a weak external electromagnetic field on the sys- 
tem under consideration. To this end, following Ref. [1], 
let us recall some main principles of the Green-function 
method. 

Consider a system that in some moment of time t is 
characterized by a statistical operator p(t). Then the 
system Hamiltonian may be written as 

H(t) =H + V(t), V(t) = J dx F 4 (x,<)^(x), 

where 7io is the Hamiltonian of an ideal atomic gas and 
V(t) is the Hamiltonian of interaction of atoms with the 
external electromagnetic field, Fi(x.,t) are the quanti- 
ties, which define the external field, and £i(x) are the 
quasilocal operators, which refer only to the studied sys- 
tem (they are not associated with the external field; see 
Ref. [J]). Note also that the summation convention over 
the repeated index i is meant. 
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Let us consider below that the influence of the external 
field on the system is weak. Due to this fact, we can 
develop perturbation theory over V(t). Thus, for the 
mean value of the arbitrary quasilocal operator a(x) in 
the linear order in V(t), one gets 

Sp/3(t)a(x) = Spwa(O) + a F (x, t), (8) 

a F (x,t) = I dt' I dx'G a ^{x-x',t-t')Fi(x',t'), 

J — oo J 

where w is the equilibrium statistical operator (the Gibbs 
operator) for the system under consideration, 



w = exp 



(9) 



Here (3 = 1/T is the reciprocal temperature in energy 
units, fi a is the chemical potential of atoms in the state a, 
and N a = / dx t)t(x)?7 q (x) is the operator of the total 
number of atoms in the state a. The temperature and 
the chemical potential are found from the following equa- 
tions: 

SpwHo — Ho, SpwN a = N a . 

The dependence of the thermodynamic potential f2 = 
fl(T,/j, a ) is defined by Spw = 1. 

The function G Q?i (x - x', t - t') in ©, 

G aU (x - x',t - t') = -i6(t - t')Spw[a{x, t),ii(x', t% 

(10) 

is the two-time retarded Green function of the opera- 
tors a(x, t), £i(x',t'), which are taken in the Heisenberg 
representation, 



a(x,i) = e lrtot a(x)e 
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The quantity 9{t) in Eq. (fT0|) is the Heaviside function. 
Taking the Fourier transforms of the quantities a F and 

« F (x,t) = j^j J dk dco e^^a F {^,oj), 
Fi(x,t) = | dk du; e- i ( fa, - kx )fi(k, W ), 
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Let us use the derived expressions for studying the 
response of a gas of alkali-metal atoms to the external 
electromagnetic field. To this end, we must replace the 
quantities Fi(x,t) in Eq. ([5]) by the vector A(x,£) and 



scalar ^(x, t) potentials of the external electromagnetic 
field. The operators a(x, t) and £, (x', t') must be replaced 
in this case by the operators of the current density j(x) 
and charge density cr(x) (see Eqs. (0])-((5|)). By the Hamil- 
tonian 7io in Eq. @ one must take the Hamiltonian TL p 
(see Eqs ([1]) and ([2])). In next calculations we use also 
the Maxwell-Lorentz equations for the average fields in 
the system (see in that case Refs. [3, 
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where the quantities cr(x, t) and j(x, t) are the mean val- 
ues of the charge and current densities (see (J3J) — ([5])), 

j(x, t) = S P/ 5(t)j(x), (t(x, t) = Sp/}(t)«r(x), (14) 

and the quantities E(x, t) and H(x, t) are the electric 
and magnetic field intensities, which can be expressed 
in terms of the potentials of an external electromagnetic 
field (see Refs. [1, @). Note that the mean values of 
the quantities <r(x, t) and j(x, t) should be calculated in 
accordance with Eq. ©. For the equilibrium state, char- 
acterized by the statistical operator ([8]), the mean values 
of these quantities vanish (see below). 

Omitting rather cumbersome calculations, we come to 
the following formulas for the mean values of the charge 
and current densities that are induced by the weak exter- 
nal electromagnetic field (see Ref. [5| for more details): 

oc 

a(x,t)= J dt' J d 3 x'(-Gi(x-x',t-t') 



x -A,(x', t') + G(x - x', t - i'Mx', , (15) 



j fc (x,t)= / dt' I d 3 x'[ -G fcI (x-x',t-0 



x -U(x', t') + G fe (x - x', t - <>(x', £')) . (16) 

The Green functions, in terms of which the expres- 
sions (Ti"5]) and (p~6|) are written, are defined by 

G(x,f) = -i6(t)Spw[a{x,t),a(0)}, 

G fe (x,t) - -*0(i)Spu;|3*(x,*),*(O)], (17) 

G fc j(x,t) = -i9(t)Spw\j k (x,t),ji(0)], 

where the charge and current density operators (see 
Eqs. (U]) and ([5])) are taken in the Heisenberg represen- 
tation. 

Let us recall that we consider the response of an ideal 
gas of hydrogenlike atoms to a perturbation by an ex- 
ternal electromagnetic field. As mentioned above, this 
consideration is reasonable in the case of a description 
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of a BEC in dilute vapors of alkali-metal atoms. Tak- 
ing into account this fact, the expressions for the Fourier 
transforms of the Green functions (see Eq. (|12[) 1 can be 
written as 



P a,P 

/ a (p-k)-/^( P ) 



£q(p) -e^(p-k) +uj + iO 



, (18a) 
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(18c) 



Here e a (p) = e a + p 2 /2m is the total energy of an atom, 
and f a (p) is the distribution function of the ideal gas of 
hydrogenlike atoms in equilibrium, which is defined by 
the relation 

Spwf] f a (p)f] a ,(p') = / a (p)<WA(p - p'), 

where 5 aa > and A(p — p') denote the Kronecker symbols. 
Bearing in mind that the unperturbed state of a gas is 
characterized by the Hamiltonian @, and using Eq. ([5]), 
it is easy to come to the explicit form for the distribution 
function f a (p), 



f a (p) = (exp{[ £Q (p)- AlQ ]/T}-l) 



(19) 



Employing the last expression and Eqs. (|4j) — (jTj) , one may 
ensure that the mean values of <r(x) and j(x) with the 
statistical operator w (see Eqs. (0) and Q) vanish. We 
recall that such mean values are not accounted in the 
derivation of Eqs. (fl"5)) and (fTB]) . Let us also note that 
the total set of the Green functions (see Eqs. (HHJ)) is given 
for a completeness of the description. Below we use only 
the expression for the scalar Green function (|18ap . 

As to Eqs. (fTS)) . it is necessary to note the following. 
Strictly speaking, these formulas are correct only in the 
case when the system is considered in a state of total 
statistical equilibrium. Thus, to produce excited atoms 
with the set of quantum numbers a in this state, one 
needs to keep matter in equilibrium with radiation. For 
that, in turn, one needs to preserve the number of pho- 
tons in the volume in which the gas of hydrogenlike atoms 
is contained (or to return the photons by the special set 
of "mirrors"). 



In the opposite case, the statistical operator ([9]) de- 
scribes a certain quasistationary state of the system, but 
not the state of total statistical equilibrium. By the term 
"quasistationary equilibrium" we mean a state of the 
system in which atoms can be found only in quantum 
states with lifetimes much longer than other character- 
istic times of system relaxation (e.g., the relaxation due 
to collision processes). For example, as these states we 
can consider the hyperfine ground states of alkali-metal 
atoms (see Ref. Q) and also states whose existence is 
stimulated by the external field (e.g., resonant laser ra- 
diation). Therefore, in the case of absence of equilibrium 
between radiation and matter, the sum over all quantum 
states a, (3 in Eqs. (fT5| must be replaced by the sum over 
these characteristic quantum states. Exactly these states 
are considered below in this paper. 

In terms of the derived Green function (I18a[) . one can 
express the permittivity of the studied system (see in this 
case also Refs. [4, 5|) as 

e~ 1 {k,u) = l + ^G{k,uj). 
This relation can be written in the form 



e(k,w) = 1 



4irk- 2 G(k,uj) 
1 + 47rfc- 2 G(k,tj) 



(20) 



Taking into account the explicit form (|18a,|l . one finds 



e 1 (k,u) = l+^-^^ f a aP (k)ap a (-k) 

P a,P 

„ /q( P -k)-^(p) 

+ i0' 



-■ ( 21 ) 

e a (p) - ep(p - k) + ' ■ 

We must note that, in accordance with Ref. Q, the for- 
mulas (fT8|) , (|20p . and (|2"Tj) can be used only in the region 
of relatively large k. Regarding the present paper, this 
criterion of applicability can be written as v\ z <C 1 (A 
is the laser wavelength and v is the density of atoms). 
In the case v}? > 1, the mentioned expressions become 
incorrect. This results from the fact that in this region 
the interaction between particles may play a significant 
role. Note that we neglect it in deriving Eqs. (fT5)) . (f^Uj) . 
and (|2ip . To take into account the interaction effects, one 
needs to derive the equations for the Green functions and 
to find a hierarchy for the correlation functions (see, e.g., 
Refs. [H Q [1] ) • In fact, the use of these procedures results 
in a redefinition of the quantity G(k,uj)k~ 2 in Eq. ([20|) 
because of the account of correlation corrections [8|, l9(. 
But, in the case when this quantity is small, Eqs. (|18|) . 
(I2TH) , and (|2~lj) and the expressions that result from them 
remain correct in the first order of perturbation theory 
even if one takes the correlations into account. Hence, 
we consider that the mentioned quantity in Eq. (|2H|) is 
small in the next calculations, 



|47rfc- 2 G(k,w)| < 1. 



(22) 



This inequality makes it possible to ignore the correlation 
corrections to the mentioned term in the first order of 
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FIG. 1: Dependence of critical temperature on the magnetic 
field intensity for the gas of alkali-metal atoms with different 
nuclear spins. The effect is caused by the Zeeman splitting 
of the lower hyperfine levels of the ground state. The gas of 
atoms is considered in the statistical equilibrium state. 



perturbation theory. We return to the discussion of this 
approximation below. 

As is known, at extremely low temperatures a Bose- 
Einstein condensate of alkali-metal atoms can be formed. 
The main peculiarity of an ideal Bose gas is that at tem- 
peratures lower than the critical temperature, T < T c , 
the chemical potential /i Q does not depend on T and 
equals the energy of the lowest level of the atomic states 
that form the BEC (see Refs. UGH), 



where by the index ao we denote the set of quantum 
numbers of the ground state. Let us note that this fact, 
in particular, results in rather strong dependence of the 
critical temperature T c on the intensity of an external 
static magnetic field (see Fig. [TJ and also Ref. [6|). 

At temperatures much lower than T c , T <C T c , we 
should consider all atoms in a BEC state. Because of 
this fact, the distribution function / Q (p) (see Eq. (ITU)) ) 
is set proportional to the Dirac delta function S(p) (see, 
e.g., Ref. 1): 

f a (p) = (2n) 3 v a 6(p), (23) 

where v a is the density of condensed atoms in the a state. 
Taking into account this fact, after integration of (T2"T|) 
over p, the formula for the permittivity of a gas in the 



BEC state (T -> 0) takes the form 

e -1 (k,a;) w 1 + -p- ^CT Q/3 (k)cr 0a (-k) 

( Va 

\UJ + As al3 - £ r + l^iap 



oj + Ae a/ 3 + s T + i-f a p J ' 

where e r = k 2 /2m is the recoil energy of an atom, and the 
quantity <r Q| a(k) is still defined by Eq. (j4|). Here we also 
introduce the natural linewidth y a p, which is related to 
the probability of a spontaneous transition between the 
a and (3 states. Note that, in our calculations, we use the 
system of units, in which the Planck constant h equals 
unity. Thus, the frequency is measured in energy units. 

As one can conclude from the relation ([24]) , in the case 
of the BEC regime, the response of the gas increases pro- 
portionally to the density of condensed atoms. This effect 
takes place due to the coherent behavior of the atoms 
in the condensate. Thus, one may say that the atoms 
respond to the perturbation cooperatively. As for the 
mathematical form, this fact corresponds to the substitu- 
tion of the wave functions of a single atom for the conden- 
sate wave functions (see Eqs. ©, and (|2"T|) - (f2~f| ). In the 
limit of zero temperatures, the mentioned phenomenon 
results in increase of collective effects, while the value 
of the corresponding linewidth (which is defined by the 
imaginary summand ijap) remains the same as for the 
isolated atom. 

It is easy to see that at frequencies close to the en- 
ergy intervals Ae a p (Ae a p = e a — ep), some peculiarities 
appear in Eq. (pH)) . Because of the small values of the 
recoil energy and the natural linewidth, these peculiar- 
ities are similar to the resonance ones. The mentioned 
behavior of the permittivity must have a strong impact 
on the dispersion characteristics of a gas and, as known, 
becomes the governing condition for slowing of the elec- 
tromagnetic pulses. 

To study the propagation of electromagnetic waves in 
a BEC we must supplement Eq. with the dispersion 
relation for free waves, which can spread in the system 

uj 2 

-^e(k,Lj)fi(k,Lj)-k 2 = 0, (25) 

where /Li(k, u>) denotes the magnetic permeability of the 
gas. This quantity may also be found in the framework 
of the Green- function formalism (see also Ref. Q). 

IV. DISPERSION CHARACTERISTICS OF A 
TWO-LEVEL SYSTEM IN THE BEC STATE 

In this section we try to simplify the description of 
the interaction of electromagnetic waves with a BEC. To 
this end, let us consider a case when the frequency of the 
external field is close to the transition frequency between 
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two given quantum states of atoms. In this case (see, 
e.g., Ref. [ll[), we can consider only resonant terms in 
Eq. (|24|) . i.e., we neglect the presence of other states of 
the energy spectrum of atoms and treat the atom as a 
two-level system. 

Let us denote by indices 1 and 2 the quantities related 
to the lower and the upper states of the two-level system 
and omit the terms that are related to other states. Then, 
Eq. (J24J) takes the form 



with 



a(k) = 47r(^i — u 2 ) 



Suj + ij 

fflff2|CTl 2 (k)| 2 



(26) 



(27) 



Here gj is the degeneracy order of the j level (j = 1,2) 
with respect to total momentum, Suj = (uj — Ae2i) is the 
detuning of the external field (\5ui\ <C Ae 2 i), and 7 = 712 
is the lincwidth related to the transition probability from 
the upper to lower state. We also neglect the so-called 
recoil energy e r (see Eq. ([24)0 because of the fact that in 
the cases studied below it is rather small in comparison 
with the characteristic quantities that give the leading 
contribution to the slowing phenomenon. 

The expression (| 2 71) can be written also in the more 
usual form (cf. Refs. @, [§]) 



c(k,w) 



1 - 



a(k) 



1 



Suj + 17 1 + C 



(28) 



where C = a(k) / (8oj + ij) . Note that the derived ex- 
pression, accurate to within a numerical factor 1/3 be- 
fore C, coincides with the Lorentz-Lorenz relation (see, 
e.g., Refs. [H, [HI]). Let us emphasize that we derive it 
in the framework of the microscopic approach. The for- 
mula d28j) results directly from Eqs. l[2T J) -([24" )) , Thus, it 
contains the microscopic parameters of the studied sys- 
tem, which are calculated in the microscopic theory. In 
accordance with Eq. (|2"3"|) . the expression (f2"B"|) carries in- 
formation about the response of the ideal gas in the pure 
Bose-condensed state. Note that in a rather dense con- 
densate there also can appear additional collective effects 
related to the interaction between atoms. Depending on 
the value of paramet er ( u\ 3 ), one should consider such 
effects (see Refs. @, S EE E3 ) ■ As it is mentioned above, 
the account of interaction effects in the region vX 3 ^> 1 
results in redefining the shift C. But let us recall that we 
consider the case of small C (see Eq. (|22p). 



\a(k)/(5u> + ij)\ < 1. 



(29) 



In some sense, this inequality results from the fact that 
we describe the system of alkali-metal atoms as an ideal 
gas. This consideration means that we neglect the in- 
teraction between atoms (see Eq. (|A1[| ). Approximate 
estimates of the applicability of this approach are given 
in the Appendix. 



To analyze the propagation properties of electromag- 
netic waves in a BEC, let us introduce the refractive in- 
dex and damping factor. We also set the magnetic per- 
meability equal to unity in the dispersion relation (|25[) . 
/i(k, u>) rs 1. It is easy to verify that this approximation 
is valid in this case. The refractive index n'(k, u>) and 
the damping factor u"(k, ui) are expressed in terms of 
the real e' and imaginary e" parts of the permittivity e, 
e = (' + ie" ', as follows: 



VVe' 2 + e" 2 + 7' 



V2 V2 
The quantities e' and e", in turn, are found from Eq. 



(30) 



Sui(Suj + a) + 7 2 
(Suj + a) 2 + i 2 ' 



7a 



(Suj + a) 2 + j 2 ' 



(31) 



It is known that the group velocity of the electromag- 
netic pulse is defined by the relation 



(32) 



n' + u>(dn' /duj) ' 



where c denotes the speed of light in vacuum. Hence, in 
the case of small energy losses and strong dispersion, one 
finds from Eq. (|3"T|) 



2c 



((<W + q) 2 +7 2 ) 2 
auj [(Suj + a) 2 — 7 2 ] 



(33) 



This expression gives the possibility of studying the de- 
pendence of the group velocity not only on the frequency 
detuning Suj and the linewidth 7, but also on the charac- 
teristic parameters of the system, which give a contribu- 
tion to the value of the parameter a (see Eq. (|27|)). 

Now we define the conditions when significant slowing 
of electromagnetic pulses in a two-level system with a 
BEC is possible without essential absorption. Consider- 
ing the limit Suj — > and employing Eq. (|3"Tj) . one gets 
the limits for the real and imaginary parts of the permit- 
tivity: 



lim e' 

Sui— »0 



7 



7 2 + a 2 



lim e" 



1 n 



7 2 + a 2 



According to these relations, the condition of small dis- 
sipation (|e"| <C e') may be written as 



7 



< 1. 



Note that this inequality defines not only the limits of 
the transparency region for electromagnetic waves, it also 
defines the upper limits of applicability of our approach 
(see Eqs. (|2"2")l and (|2"9"|)). Really, as it is easy to see, in 
the case \a\ ~ 7, the second term in the denominator 
of Eq. (f28|) (the so-called Lorentz-Lorenz shift C) can 
become of the order of unity. But, as mentioned above, 
we assume that this term is small (|C| <C 1). 
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We need also to use the condition defining sufficient 
pulse slowing (the strong dispersion condition), v g <C c, 
which, in accordance with ((33|) . take the form 



Aggijaj 
2 7 2 



» 1. 



(34) 



In a system with defined energy structure of atoms 
(Ae, 7 = const) the only parameter that may be var- 
ied is the level occupation difference (y\ — v%), which is 
included in the parameter a (see Eq. (|27)), Taking into 
account this fact and the above inequalities, one gets 



7 M 
-j— — < — < 1. 

Ae 2 i 7 



(35) 



This inequality characterizes the region of densities (pa- 
rameter a) in which the phenomenon of significant slow- 
ing of electromagnetic waves in a BEC may be observed. 

Next, on the basis of the derived formulas, we consider 
two cases of slowing of weak electromagnetic waves in 
a BEC. The first one is concerned with electromagnetic 
pulses in the optical region of the spectrum (resonance 
D lines of alkali- metal atoms), and the second one is re- 
lated to microwaves (hyperfine transitions of the ground 
states of alkali metals). 



V. SLOWING OF ELECTROMAGNETIC 
WAVES IN THE OPTICAL RANGE IN A BEC 

As was mentioned above, a significant slowing of elec- 
tromagnetic pulses was observed in a condensed gas of 
sodium atoms in Ref. T]. In this experiment, for the 
observation of ultraslow light pulses the effect of electro- 
magnetically induced transparency (EIT; see Ref. [lj| for 
details) was used. For realization of this effect one needs 
to use an additional laser of large intensity besides the 
laser pulse, which is slowed in the experiment. But we 
should note that the description of this experiment in the 
framework of the Green-function formalism is impossible. 
The main reason is that for EIT support it is necessary 
to use a high-power (coupling) laser. Thus, the influence 
of such a laser on a BEC cannot be described correctly 
within the linear response theory. But let us emphasize 
that the use of a high-power laser is not a necessary con- 
dition for the existence of the slowing phenomenon. 

Therefore, we consider the above microscopic theory 
in the case of slowing of weak light pulses in a sodium 
BEC at temperatures T <C T c in the absence of EIT. In 
particular, in the framework of the developed approach 
we study the slowing phenomenon, when the laser fre- 
quency is close to the transition from the ground state 
(I = 0, 35*1/2 state, F — 2 for sodium) to the excited 
one (I = 1, 3P 3 / 2 state, F' — 1). Note, that the dipolc 
transition 1 — > 2 between these states is allowed. 

First of all, we emphasize the following. The matrix 
element of the charge density U\ 2 (see Eq. @) contains 
the wave functions of a hydrogenlike atom <y9i. 2 (y), which 



have a sharp maximum at the point y = 0. Due to this 
fact we can expand the exponents exp [z(m p /m)ky] and 
exp [— i(m e /m)ky] into series of perturbation theory over 
(ky) <C 1. Then, in the first order of perturbation theory, 
one gets 

ojaHk) w ikdi 2 , 

where di2 is the atomic dipole moment, which is related 
to the dipole transition 1 — ► 2: 
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eri2, ri2 = J dy <p*(y)<f2(y)y- 



(36) 



For the quantitative estimates we take for the sodium 
atom |cr 12 (k)| 2 w k 2 d 2 /3 that results in (see Eq. (j27| ) 
a = (4tt/3)Vi - v 2 )d 2 . Taking d 2 » (3.52er ) 2 SW' [11, 
where ro is the Bohr radius and Sff> is the relative 
strength of the F — ► F' transition, Ae2i ~ 2.1 eV, and 
the linewidth 7 = SpF^e/^ (r e = 4.1 x 1(T 8 eV is 
the natural linewidth [lq : below we consider F = 2 — > 
F' = 2 transition, S22 = 1/4). Since we consider the 
propagation of weak electromagnetic waves through the 
condensate, we assume that the density of atoms in the 
excited state is negligibly small. In other words, we treat 
the occupation of the excited level is small due to low in- 
tensity of the propagating pulse (y 2 <C v\)- Next, using 
Eq. (|35p . one can obtain the following inequality, which 
defines the conditions of slowing in condensed sodium 
vapor at T <C T c : 



3 x 10 5 < v < 10 14 cm" 3 . 



(37) 



From this inequality it is easy to see that the density 
value of condensed atoms in the case of real experiments 
with BECs {v = 5 x 10 12 cm" 3 ; see Ref. Q) lies in the 
limits of the derived interval. 

Note that the main parameter that can be mea- 
sured in experimental conditions is the relative inten- 
sity of the transmitted light. This quantity for a ho- 
mogeneous medium is defined by the expression I / Iq = 
exp [— e"(u})kL], where L is the transversal size of the 
atomic sample. Taking the density of condensed atoms 
in the ground state v = 5 x 10 11 cm~ 3 , k = 1.066 cm -1 , 
L = 0.004 cm, one finds I/I w e~ 2 - 2 w 0.11. This 
allows to state that in the limit of zero temperatures, 
T <C T c , ultraslow light pulses of weak intensity may 
be directly observed for these parameters (according to 
Eq. (f3~4"|) . Vg/c « 10~ 6 ). However, a significant part (up 
to 89%) of the pulse is absorbed by the medium. 

It is easy to see that the intensity of the transmitted 
light has a rather strong dependence on the system pa- 
rameters, which make a contribution to the exponent. 
Thus, the situation with the slowing phenomenon may 
completely change, when temperature effects are taken 
into account. These effects may increase the value of 
e"(u>) resulting in the strong reduction of the relative 
intensity of the transmitted light. For example, with in- 
crease of e"(uj) by several times the intensity reduces to 
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FIG. 2: Scheme of energy levels of a sodium atom. For the 
sake of simplicity the levels are drawn not to scale, and the 
hyperfine structure of the excited state is not shown. 



less than l%. It is possible to take into account the influ- 
ence of temperature effects on the propagation properties 
in the framework of the Green-function formalism. How- 
ever, this is a separate problem, which is not considered 
in the present paper. 

Note also that in the case described above the group 
velocity of the pulse is negative (anomalous dispersion in 
the case of a normal population, v\ > v-i). We return to 
the discussion of this fact in Sec. I VII 

Next, let us consider light slowing in a multicompo- 
nent BEC employing the developed theory. Nowadays, 
the realization of a multicomponent Bose condensate in 
experimental conditions is rather achievable (see, e.g., 
Ref. The formulation of this problem is interesting 

for the following reasons. As is known, in an external 
magnetic field the Zeeman splitting of levels occurs, and 
the spacing between different components of the multi- 
plet Ae mag is proportional to the magnetic field intensity. 
Thus, in a weak external magnetic field this energy dif- 
ference is much smaller than the difference between the 
nonsplit levels of the hyperfine transition Ae n f- If the 
atoms can occupy different magnetic states of the mul- 
tiplet, we can talk about the existence of a multicompo- 
nent BEC. In such a condensate the effects related to the 
resonant absorption for some defined frequencies may be 
significantly smaller than in the case described above. 

As an example, we consider a two-component conden- 
sate of sodium atoms in the quantum states 3Si/2, F = 2, 
mp = —2 and 351/2, F — 2, rap = —1 (below we de- 
note them with subscripts 1 and 2, respectively). Now, 
let us study the response of the system to laser radia- 
tion, which is detuned relative to the transition between 
these two states and the state 3P3/2 (denoted below by 
the subscript 3), as shown in Fig. [2] 

By using Eq. , we study the dependence of the dis- 
persion characteristics of an ideal gas of sodium atoms 
in a BEC state on the frequency of external radiation. 
For simplicity, let us treat as equal the transition prob- 
abilities from the upper state to one of the lower states, 
731 = 732- We also consider the densities of atoms in 
two lower states v\ = v% = 5 x 10 12 cm~ 3 and 013 = 023. 




0.0 U 1 L 
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c ' mag 

FIG. 3: Dependencies of the refractive index, damping factor, 
and transmitted light intensity for the two-component Bose 
condensate. 



In fact, these assumptions are made only for convenience, 
but they are not necessary in principle for numerical eval- 
uations. Let the energy of splitting be Ae mag = 8731. 
Taking all other parameters the same as in the above 
estimates, one finds the dependencies of the dispersion 
characteristics, as illustrated in Fig. [3] From these graphs 
we conclude that in such a system the signal may slow 
down approximately to 264 m/s. 

At first sight, the last result contradicts a well-known 
statement from classical optics (see, e.g., Ref. 18]) that 
in a two-level gas a large index of refraction is usually 
related to a high absorption ratio. Hence, on atomic res- 
onances the signal is mostly absorbed. However, this con- 
tradiction is only apparent. Let us emphasize one more 
time that in the case described above the real part of the 
susceptibility x' = ( e ' — l)j which gives the main con- 
tribution in the change of the refractive index, is rather 
small in magnitude (as one can conclude from Fig. [3J 
Ix'l i$ 1CP 2 ). The imaginary part x" = e " \ which defines 
the absorption ratio, is rather small too, x" ^ 10~ 2 (in 
the central part x" ~ 10~ 3 ). In other words, in the con- 
sidered case, the matter is not dense for laser radiation 
(the refractive index is close to unity), but has a strong 
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dispersion due to closely situated resonances. This is the 
reason why the group velocity of the pulse may be suffi- 
ciently reduced, while its absorption ratio remains rather 
low. 

Let us note that in the regions of anomalous disper- 
sion in Fig. [3] (where the derivative of the refractive in- 
dex is negative) there is good agreement with the results 
obtained in Ref. [12], in which the authors studied the 
dependence of the refractive index on the laser intensity. 
The agreement takes place in the limit of low light inten- 
sity. In fact, we consider only this limiting case in the 
present paper. 

As for the obtained dependencies as whole, they are 
similar to the dependencies that are observed in the EIT 
regime [Hj]. One can, by analogy with Ref. [l|, state 
that the phenomenon of simultaneous manifestation of 
low group velocity value and low absorption ratio in the 
central part of Fig. [3] is the result of destructive interfer- 
ence between two absorption-emission paths, which cor- 
respond to the transitions 1 — > 3 and 2^3, respectively. 
The difference is that the spacing between these transi- 
tions equals not the Rabi frequency value Qr, but the 
energy of splitting of two occupied magnetic sublevels 
A£ mag . If we compare the dependence of light absorp- 
tion (the lowest gra ph i n Fig. [3]) with the analogous one 
in the EIT scheme [15| . we can see that the intensities 
of the transmitted light in the central part are nearly 
the same (I/Iq « 0.3). However, in the EIT case, the 
maximum is sharper due to the nonlinear nature of the 
effect. 

Let us emphasize that in our case the energy spac- 
ing between magnetic states must be larger than the 
linewidth of the transitions. In other words, these sub- 
levels cannot be situated arbitrarily near each other. 
Therefore, the intensity of the transmitted light depends 
directly on the energy spacing Ae mag and on the density 
of the condensed states in the condensate v\^. At the 
same time, the group velocity of the pulse also strongly 
depends on both these characteristics. Thus, if one tries 
to change the characteristics of the system for much 
stronger reduction of the group velocity value (increase 
in the density of atoms or decrease in the energy spac- 
ing Ae mag ), the ratio of absorption will correspondingly 
increase (see also Ref. [l8|]). In our case, the absorption 
ratio is the main limiting factor for the achievement of 
lower values of the group velocity. 

As for the effect of electromagnetically induced trans- 
parency, in the experiment the Rabi frequency (which is 
analogous to the energy Ae mag in our case) is usually set 
to be of the same order as the linewidth. Thus, the ef- 
fect of slowing due to the Zeeman splitting is somewhat 
weaker (v g « 264 m/s in comparison with v g w 32 m/s in 
Ref. [l[). However, the authors hope that the observation 
of this kind of slowing is rather possible in the present 
experimental conditions. 



VI. SLOWING OF MICROWAVES BY 
HYPERFINE GROUND STATES OF 
ALKALI-METAL ATOMS 

As mentioned above, the proposed microscopic the- 
ory is correct when the occupation of the excited atomic 
states is stimulated by an external electromagnetic field 
(the case that is studied in Sec.|V|, or when the lifetimes 
of the excited states are much larger than the relaxation 
times. Note that for alkali-metal atoms all sublevels of 
the hyperfine ground state can be treated as such long- 
living states. This follows from the fact that the dipole 
transitions are forbidden for these states. Waves whose 
frequency corresponds to the energy difference between 
the levels of hyperfine structure are usually called mi- 
crowaves. Since there are many trapping-related experi- 
ments with ultracold atoms prepared in different hyper- 
fine states, the problem of the possibility of slowing mi- 
crowaves becomes rather significant. 

By considering the example of cesium atoms in the 
BEC state, we try to find the conditions in which, in 
our opinion, the slowing phenomenon may be observed 
experimentally. To this end, according to the developed 
approach, let us show that the region defined by Eq. ([35]) 
exists. Note that a brief description of the application 
of the microscopic approach for the description of mi- 
crowave slowing may be found in Ref. [191 ]. 

It is known that alkali-metals ( 133 Cs, in particular) 
do not have a dipole moment in the ground state. This 
agrees with Eq. (f3"6"|) . from which one can find that the 
dipole moment d equals zero for the wave functions of 
the ground state. Thus, to describe the slowing of mi- 
crowaves on the levels of hyperfine structure with for- 
bidden dipole transitions, we should expand the charge 
density matrix element o\i (see Eq. ([J}) to the second 
order in (ky) <1. As a result, one gets 

~(*ro) 2 , (38) 

where r is the atomic radius (for the cesium ground 
state tq ~ 2.6 x 10~ 8 cm [20(). Let us consider that the 
levels are degenerate, g x = 7, g 2 — 9 (gj = 2Fj + 1, 
where Fj is the total spin of an atom in the j hy- 
perfine state, j = 1,2), and also set k = (Ag2i/c) 
with Ae2i ~ 3.8 x 10~ 5 eV (microwaves with frequency 
9.2 GHz). We take the natural linewidth 7 of the upper 
hyperfine level equal to the value that correspond to the 
anticipated accuracy (10~ 16 ) in "cesium fountain clock" 
experiments [2l[, 7 w 3.8 x 1CP 21 cV. Then, on the basis 
of Eqs. ([2"7]) and ([55]) , one may get the following region 
where the slowing phenomenon may be observed: 

1CT 1 < \y t - y 2 \ < 10 15 cuT 3 . (39) 

It is easy to see that the effect becomes greater when the 
density difference increases. When it reaches the upper 
limit of (|39p . damping effects prevail in the system and 
the signal is mostly absorbed by the medium. We stress 
that the region of densities ([3"9"[) . in our opinion, looks 
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FIG. 4: Level scheme of cesium atom in an external mag- 
netic field. The first sign in brackets corresponds to the total 
spin F, the second one to its projection rriF- For the sake of 
visibility, the figure is not drawn to scale and only the most 
relevant levels are shown. 



convenient from the standpoint of experiments with the 
BEC regime in cesium (see, e.g., Ref. [22j], in which ce- 
sium atoms with the density v = 1.3 x 10 13 cm -3 were 
used). 

Now we note the following. As it is easy to see from 
Eq. (|33|) . in the limit 8lo — > 0, the sign of the group 
velocity v g depends directly on the sign of the quan- 
tity a, which, in turn, depends on the sign of the dif- 
ference (vi — V2). In other words, it depends on whether 
the population is normal or inverse. In the case of nor- 
mal population (y\ > 1/2), the group velocity of the signal 
is negative. It is traditionally considered that the group 
velocity for the transparent matter is positive, and it can 
be negative only in regions in which the pulse is mostly 
absorbed. But, from our theory one can conclude that, 
according to Eq. ([33]) . the signal can propagate in the 
system with rather small dissipation (i.e., the matter is 
transparent) and rather slow velocity. Note that the ob- 
servation of electromagnetic pulses with negative group 
velocity is not so abnormal. The existence of this kind 
of phenomenon for physical systems when the wave fre- 
quency is close to atomic (or molecular) resonances was 
pointed out and studied in many works (both theoreti- 
cal [H, HI| and experimental [Ifl [26[ ) . In the case of in- 
verse population (y\ < V2), more "usual" situation takes 
place because the group velocity of the slowed pulse is 
positive. 

So far, we assumed that the external magnetic field is 
absent, i.e., we treated the hyperfine levels as degener- 
ate in the quantum number rap. However, an external 
magnetic field (such situation occurs in most experiments 
with BECs) results in a more complicated picture. Each 
hyperfine level splits to additional sublevels with differ- 
ent total spin projections mp (e.g., for the cesium atoms, 
see Fig. [4]). As is known, in a weak (and also in a strong) 
external magnetic field the energy spacing between such 
sublevels is proportional to the magnetic field intensity 
(the Zeeman and Paschen-Back effects, respectively). 



Let us assume that the magnetic sublevels of the up- 
per level of hyperfine structure are kept away from each 
other (the linewidth of the transition is much smaller 
than the energy difference between sublevels Ae mag , as 
shown in Fig. |4|). Then, a signal tuned exactly to the 
transition between the occupied lower state and one of 
the upper states (narrow black arrowed lines in Fig. 3|) 
may be slowed down, as described above. The only dif- 
ference from the previous situation is that the states are 
no longer degenerate (51 = 52 = 1) and the linewidth jj 
in Eq. (|2l)|) corresponds to the transition probability from 
the chosen sublevel with the projection j = (—F\ F') 
to the occupied lower state. Taking all other values the 
same when deriving the inequality (|39[) , one can find the 
inequality in the presence of the external magnetic field 



10 < \ui - v 2 



< 10 17 cm- 3 . 



(40) 



In the experiment [22j . condensed cesium atoms with a 
peak density of 1.3 x 10 13 cm -3 were kept in a trap. Thus, 



according to Eq. (|40|) . one finds that in such a system 
ultraslow microwaves may be observed. Moreover, direct 
calculations show that the group velocity of the signal 
(see Eq. (1331) ) propagating in the mentioned system is 
close to 3 x 10~ 4 cm/s. 

Note that here we must make important remarks re- 
lating to the possibility of experimental observation of 
ultraslow microwaves in a BEC. In fact, it is impossi- 
ble to make a direct measurement of the group velocity 
value in the conditions of a real experiment. It is usu- 
ally calculated by measuring the time of pulse delay. The 
characteristic dimensions of the pulse (and, specifically, 
the wavelength) must be of the order of the atomic sam- 
ple size, in order for it to propagate in the system. In 
the present experimental conditions this relation is not 
always satisfied. Usually, the characteristic sizes of a con- 
densate are of the order of a millimeter (or even submil- 
limeter), whereas the wavelength, related to the hyperfine 
transition, is of the order of centimeters (e.g., for cesium 
in the absence of an external field Ahf = 3.3 cm). 

In view of the aforesaid, there are two possible ways 
to avoid this problem. The first one, which is rather 
difficult at this moment, is to use other configurations of 
magneto-optical traps (and, probably, other cooling tech- 
niques) that may allow production of condensed clouds 
with centimeter sizes. The second one, and probably 
more realistic, is to use a rather strong bias field, which 
results (according to the Paschcn-Back effect) in increase 
of the spacing between levels and, consequently, in the re- 
duction of the resonance wavelength. At the present mo- 
ment, such fields are completely accessible, although they 
must be rather strong. For example, to reduce the clock 
transition wavelength to Ahf — 2 mm in cesium, the field 
intensity must be of the order of 45 kG, which is quite 
achievable from the standpoint of the trapping-related 
experiments (see Ref. [27| )• Let us note also that, corre- 
sponding to the results p (see also Fig.[T]), the presence 
of the additional magnetic field may sufficiently raise the 
value of the critical temperature. 
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Note that, in order to observe the effect for the limiting 
case considered while deriving Eq. (|55|) . one needs to use 
a microwave signal with detuning duj s much smaller than 
the level lincwidth: 



Suj s < 7j . 



(41) 



In a real experiment it is rather difficult to to satisfy 
such a condition. This is due to the fact that the natural 
linewidth jj of the hyperfine structure levels is rather 
small (dipole transitions are forbidden and the transitions 
come from the higher order effects). 

In this case, the detuning of a maser between close- 
lying magnetic sublevels of the hyperfine ground state 
structure may save the situation. At the same time, we 
must satisfy the condition 

where Ae mag is the energy difference between these sub- 
levels. For example, a microwave signal that is detuned 
relative to the transitions between two neighboring states 
of the upper hyperfine multiplet and the occupied lower 
state (wide gray arrowed line in Fig. [4|) may satisfy this 
condition at certain intensities of the magnetic fields. 

Let all atoms be condensed in the lower (occupied) 
state with the density v and energy e - Then, it is not 
difficult to come to the following expression for the per- 
mittivity (see Eqs. (HD and l[27|)): 



F' 

• £ 

J = -F> 



1 



(42) 



where Aej — (e\F',m F t) ~ £ o)- The derived expression 
defines the dependence of the refractive index on the fre- 
quency of the microwave radiation. Therefore, it also 
defines the conditions for pulse slowing (see Eqs. ([30)) 
and (|32[) ). It is easy to see that in this case the condition 
of wave slowing (cf. Eq. ([55]) takes the form 



Ae„ 



As 



4a 



hf 



Ae 



Ae r 



mag 



20 



(43) 



The structure of the 133 Cs upper hyperfine state in the 
presence of magnetic field is schematically shown (not 
to scale) in Fig. [5] From this graph, one may conclude 
that the energy difference Ae mag is small either in the 
region of weak magnetic fields or in the region of levels 
intersection (strong magnetic fields). 

In the presence of an external magnetic field the de- 
pendence of the refractive index has a rather complicated 
form due to the large number of magnetic states (nine 
sublevels of the upper F = 4 state; see Fig. Hence, it 
is more convenient to show such dependence in a figure. 

The dependencies of the refractive index and damping 
factor on the frequency are shown in Fig. [6l Notice that 
in the central part [5uj = 0, corresponding to the wide 
gray arrowed line in Fig. U) the refractive index has a 
slope with steepness dependent on the interval Ae mag = 
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FIG. 5: Energy spectrum of the upper (_F=4) hyperfine 
multiplet of the 133 Cs ground state. The region of inter- 
section of the two upper levels corresponds to the inten- 
sity B CI « 583 kG. 



e |4,4) — £|4,3)i which, in turn, depends on the magnetic 
field intensity. This means that a microwave signal that is 
detuned in this way may be sufficiently slowed. Moreover, 
its group velocity depends directly on the magnetic field 
intensity. 

The direct calculations show that, e.g., for an ideal 
gas of condensed cesium atoms (v = 1.3 x 10 13 cm~ 3 , 
7j = 3.8 x 10~ 21 eV) with the energy of splitting Ae mag = 
5.7 x 10~ 20 eV (Ae mag /7j = 15, such as shown in Fig.JHJ), 
the pulse can be slowed down to 0.01 cm/s. But we em- 
phasize that in this case (weak bias field) the pulses with 
wavelengths of the order of centimeters may be slowed 
in a BEC. Thus, one must use condensed clouds of the 
same size, which is, as mentioned above, hardly achiev- 
able from the standpoint of present experimental condi- 
tions. 

In the case of strong magnetic fields (close to the region 
of intersection of levels of the Paschen-Back splitting; see 
Fig. [5]) , the graphs of the refractive index and damping 
factor have qualitatively the same dependence, as shown 
in Fig. [51 Moreover, the wavelength of the corresponding 
radiation in these conditions is of the order of submil- 
lfmetcrs, which is quite acceptable from the standpoint 
of a real experiment. But, as follows from the estimates 
that are based on Eqs. ([4*2")) and ([4*3")) , the conditions for 
slowing of electromagnetic waves are satisfied only in a 
region of magnetic field intensity B cr that approximately 
equals 583 kG. In this region of bias field one can ef- 
fectively control the ultraslow microwave phenomenon in 
a BEC. Present scientific papers note the achievement 
of ultrastrong intensities of magnetic fields (e.g., of the 
order of 600 kG; see in that case Ref. [28[ ), but the au- 
thors of the present work, being theorists, cannot state 
whether it is possible to use such fields in BEC-related 
experiments. 
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FIG. 6: Refractive index and damping factor dependencies 
for cesium atoms in the BEC state in the region of frequen- 
cies close to the transitions between magnetic sublevels of 
hyperfine splitting. Left and right peaks correspond to the 
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tively. The parameter used for the Zeeman splitting energy 

A£ mag /7.J =15- 



VII. CONCLUSION 

Thus, we propose a microscopic approach for the de- 
scription of the linear response of an ideal gas of hy- 
drogenlike atoms in the BEC state to perturbation by 
an external electromagnetic field. The use of this ap- 
proach allows us to find the permittivity of dilute gases 
of alkali-metal atoms in a BEC state as a function of 
the microscopic characteristics of the system at temper- 
atures T < T c . 

The main advantage of the considered approach is that 
it comes from a microscopic description, which is based 
on the approximate formulation of the second quanti- 
zation method for systems that consist of bound states 
of particles [3| . Another advantage is its convenience for 
studying the dispersion characteristics of gases regardless 
of the type and number of quantum states in the system. 
As a possible drawback one can refer the fact that in 
the form in which it is formulated in this work it cannot 
be correctly used for the description of the response to 
a strong perturbation (as may probably be done in the 
approaches [12|, ll4[) owing to the use of linear theory. In 
other words, this approach cannot be used without con- 
siderable modification to systems in which strong laser 
radiation is used (e.g., the coupling field for the provision 
of EIT) . Note that we neglected also the interaction be- 
tween atoms (in view of the consideration of rather dilute 



gases) and the influence of temperature effects. Temper- 
ature effects in the theory developed may be taken into 
account, but this represents a separate problem. 

The developed approach is demonstrated on the ex- 
ample of the reduction of the speed of light in a BEC of 
sodium vapor with parameters close to the experimental 
ones [l|. It is shown that in some conditions the signal 
may be sufficiently slowed in the absence of an additional 
strong (coupling) field. A strong field is used only to pro- 
vide electromagnetically induced transparency leading to 
the reduction of the damping factor value. Thus, it is 
not a necessary condition for the observation of ultra- 
slow pulses. 

The following cases are considered: the optical signal 
is tuned to the resonance frequency and the signal is de- 
tuned relative to close transitions in the two-component 
condensate. It is shown that in these situations ultra- 
slow pulses may be observed. The conditions for slow- 
ing, which depend on the macroscopic characteristics of 
the system, are found. The dispersion dependencies are 
studied also. 

The approach is applied for studying the propagation 
of microwaves in a BEC. It is shown that the microwave 
signal, which is tuned to the transition between different 
hyperfine levels of the ground state, can propagate in 
the system with rather small energy loss. It is pointed 
out that the sign of the group velocity of the propagating 
pulse depends directly on the occupation difference of the 
states. 

The influence of an external static and homogeneous 
magnetic field (taking into account the Zeeman and 
Paschen-Back effects) on the slowing parameters is stud- 
ied. Considering cesium atoms, it is shown that the pulse 
can be slowed down to extremely small values in some 
conditions. 



APPENDIX A: ON THE CORRECTNESS OF 
THE IDEAL GAS APPROXIMATION 

Let us now say some words relating to the limits of 
applicability of the developed approach. Note that we 
consider the condensed gas of alkali-metal atoms as ideal. 
In other words, we assume that the interaction between 
atoms is small in comparison with the binding energy 
and in comparison with the interaction of particles with 
the external field. 

To estimate the influence of effects related to the inter- 
atomic interaction, let us consider the Hamiltonian 7ii n t , 
which is introduced in Eq. {!]). As mentioned above, the 
microscopic expressions for it are found in Ref. Q in the 
framework of the developed formulation of the second- 
quantization method. It is also shown that in the low- 
energy approximation (or in the case of low tempera- 
tures) this operator can be written as 

Him = / dx dy v a0 ^ s (x-y)f}^(x)f}l(y)fj J (y)fjs(x), 



(Al) 
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where the Einstein summation convention over indices 
a,/3,7,(5 is meant, and the kernel v a p.-ys( z ) has the form 

v a /3; 7 s{z) = -g [^(dasdfa) - 3(zd Q 5)(zd (37 )] , (A2) 

in which the quantities d a p are the matrix elements of the 
dipole moment operator of a single atom. Such matrix 
elements can be expressed in terms of the atomic wave 
functions <p a (y) in different quantum states a,/3, 



where 



la/3 



(A3) 



Thus, one can see that in this case the Hamiltonian (|A1[) 
corresponds to the dipole-dipolc interaction of atoms. 

It is easy to verify that in the case of the interaction 
of atoms in the ground state the expression (|A1|) is inap- 
plicable. It results from the fact that the dipole moment 
of the hydrogenlike atom in the ground state is equal 
to zero. One may find the Hamiltonian of interaction 

Ti-lnt °f the atoms in the ground state from the second- 
order terms of the perturbation theory over the interac- 
tion (jAip . If we assume that this interaction is small 
in comparison with the binding energy |J2J|, we get (see 
Ref. 1) 

^int =\J d*dy ^ )(x-y)r7t o (x)ryt o ( y ) 7 ) ao ( y ) 7 ) Qo (x), 

(A4) 

where ry a (x) and f] ao (x) are the creation and annihi- 
lation operators of hydrogenlike (alkali-metal) atoms in 
the ground state with the set of quantum numbers ao 
at the point x, respectively. The quantity i/°)(z) can be 
expressed by the relation 



E'l W /3A;a a ( Z )| 2 

» lJ W = , 



(A5) 



where v/3\^ aoao (z) is defined by Eq. (|A2|) . and e ao is the 
ionization energy of an atom. The prime in the sum 
means that the terms with (3 = and A = «o are omit- 
ted. Hence, one can conclude that the Hamiltonian (|A4|) 
corresponds to the well-known Hamiltonian of the van 
der Waals interaction (see, e.g., Ref. [29|). 

We must note that it is convenient to replace the op- 
erators (jAll) and (IA4[) by effective Hamiltonians (see in 
this case Refs. [7|, |30(). Thus, the effective Hamiltonian 
of the dipole-dipole interaction takes the form 

^inf 5 = U af 3- n s / dx ?) a (x)^(x)?7 7 (x)i7 < 5(x) (A6) 



with (see also Ref. [3l|) 



4tt 



— — dagdpy. 



(A7) 



In the case of interaction of atoms in the ground state, 
the Hamiltonian (|A4j) can be replaced by 



^int = \ u o J dx ?7 ao (x)4 (x)^o( x )»?«o( x )) ( A8 ) 



U 



Anh a s 



(A9) 



and a s is the s-wave scattering length (see, e.g., Ref. [30|). 
Note that the quantities U a p- tl s and Uq may be defined by 
Eqs. (|A7[) and (|A9[) only in the case when the quantum 
fluctuations, which are related to the elementary excita- 
tions in the system, are negligibly small [3(3]. 

To find the approximate criteria of applicability of our 
theory, let us consider a dilute weakly nonideal Bose 
gas, which interacts with a weak external electromag- 
netic field. It is known that the standard thermodynamic 
perturbation theory is inapplicable to such a gas in the 
region below the critical point T c . It results from the fact 
that there appear diverging terms in the region of small 
momenta (see in this case, e.g., Refs. 0, H(|)- Therefore, 
one needs to use a special perturbation theory. Such 
a theory was constructed by Bogoliubov in [32j ■ Within 
this theory it is shown (and rigorously proven in Ref. (33|) 
that the creation and annihilation operators may be con- 
sidered as c numbers, 



(x) 



Vao ( x ) 



where <fi is some phase, which is irrelevant in our con- 
sideration. To account for the existence of the overcon- 
densate particles (the density of which is much smaller 
than the density of the condensate particles) one needs 
to introduce a gas of elementary excitations, or quasipar- 
ticles with an energy spectrum, which is now called the 
Bogoliubov spectrum. The special perturbation theory, 
which was proposed by Bogoliubov, removes the diver- 
gences of the correlation functions in the region of small 
wave vectors [3(| ■ In addition to the mentioned spec- 
trum, this theory allows one also to find the ground state 
energy of the weakly nonideal Bose gas, its chemical po- 
tential (which has a nonzero value in case of interaction) , 
and the corrections to these values, which result from 
the existence of quantum fluctuations. For example, the 
contribution of interaction effects to the ground state en- 
ergy £ can be written as (see Ref. [3T 



£o = U 



2V 



1 



128 



{vat) 



3\l/2 



(A10) 



where U is the physical coupling constant (see Eqs. (|A7[) 
and (|A9|) ). N is the overall number of atoms in the weakly 
nonideal Bose gas, V is the system volume, v is the 
density of particles, and a s is the scattering length (see 
Eq. (|A9[0 . The second summand in brackets in Eq. (|A10[) 
characterizes the quantum fluctuation contribution to the 
ground state energy of the gas. The formula (|A10[) is cor- 
rect, when the so-called gas parameter ^|a s | 3 is small, 



< 1. 



(All) 



Below we assume that the relation (|A11[) is satisfied. 
Therefore, in the following considerations we do not take 
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into account the contribution of the quantum correla- 
tions to the characteristics of the weakly nonideal Bosc 
gas interacting with the external radiation. 

We note that the formalism developed in the present 
paper can be applicable to the system under considera- 
tion only in the case when the interatomic interaction is 
small compared to the interaction with the external field. 
Thus, let us compare the energy related to the interaction 
with the weak external radiation and the contribution of 
the interaction between atoms to the ground state en- 
ergy £ . If we neglect the existence of the quantum fluc- 
tuations in the system, this contribution, in accordance 
with Eq. (|A10[) . takes the form 



(A12) 



where U a is the coupling constant, which characterizes 
the effective Hamiltonian of interaction between atoms 
with the set of quantum numbers a (see Eqs. (|A6I) — (|A9|) ) . 
and u a is the density of such atoms in the condensate (by 
the index ao we denote the set of quantum numbers of 
the ground state). The energy related to the interaction 
with the external electromagnetic field is characterized by 
the Hamiltonian V(t) (see Eq. (J3j) ) . Its mean value £e 
can be expressed in the framework of the Green-function 
formalism (see, e.g., Refs. [1, H[). But, one can estimate 
it also in a different way. This energy must be of the 
same order of magnitude as the energy, which is pumped 
into the system by the external field [J]. As mentioned 
above, the expression for the energy that is absorbed by 
the system from the external field can be written as (see 
Eq. (JUD) 



£e 



8^' 



(A13) 



where E is the amplitude of the electric field in the pulse, 
£ = {1 — exp [— e"(uj)kL]} is the part of the energy that 
is absorbed by the system in the case of propagation of a 
single pulse, and L is the characteristic linear size of the 
atomic cloud. 

The ideal gas approximation is correct when the energy 
£q is small in comparison with £e, 



^ U a v a v ao < — f . 

07T 



(A14) 



In Sec. IIVI we study the dispersion characteristics of a 
two-level system in the BEC state. There we consider 
that the density of atoms in the excited state (denoted 
by the set of quantum numbers a>i) is much less than 
the density of atoms in the ground state (denoted by the 
set of quantum numbers ao), v ax -C v ao - The number 
of atoms in the excited states can be estimated with the 
assumption that all energy (IA13P that is absorbed by the 
system is spent on the excitation of certain states a±. If 
the intensity of the electromagnetic field is rather small, 
one gets (see, e.g., Ref. 34]) 



£2 



(A15) 



From the inequality v ai <C v ao we have 
E 2 



(A16) 



Therefore, using Eq. (|A15|1 we come to the following re- 
lation for the contribution of the interaction (|A12[) in the 
ground state energy (see also Eqs. (|A7|) . (|A9[) ): 



i 2 a s 2 47T 2 E 2 



Thus, the inequality (|A14[) can be written as 



Anh 2 



47r . |2 E\ E\ 



(A17) 



This relation defines the limits for the intensity of the ex- 
ternal electromagnetic field and the density of condensate 
particles in the ground state. 

To find numerical estimates, we take the value of the 
field intensity from the experiment [l[, in which a laser 
with power density 5 mW/cm 2 and A =589 nm was used. 
The scattering length a s for a number of alkali-metals in 
the ground state is of the order of a few nanometers [3(| • 
Taking the density v ao 
easily obtains 



ATrh 2 a : 



8 2 



47T 2 E?_ 



E^ 

8vr 



10 cm 6 and £ « 0.5, one 

- 10~ 12 erg cm -3 , 
10~ 14 erg cm~ 3 , 
10~ 6 erg cm -3 . 



It is easy to see that in this case the inequality (|A17|) is 
correct. Let us note that the relation (|A16[) is satisfied 
also with a rather large "reserve" . 

Now we can conclude that the relations ([22]) and ([29]) 
(see also Eqs. ([37]) and (|30"]0 and the inequalities (|A14|) . 
(|A16|) , and (|Al7j) define the limits of applicability of our 
theory for the description of the response of a condensed 
gas of alkali-metals to a perturbation by a weak external 
field. For completeness of the description, one should 
also add the condition 

E 2 

which characterizes the smallness of the energy related 
to the interaction with the external field in comparison 
with the characteristic energy that may be enough to 
ionize atoms in the condensate. 

Note also that here we do not consider the interac- 
tion of closely located radiating dipoles. This can be 
explained by the fact that the effects produced by this 
interaction are quadratic in the density of atoms in the 
excited state (dipoles), and the density of dipoles in di- 
lute vapors of alkali-metals is low (the average distance 
between dipoles is much larger than the characteristic 
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wavelength). First, the experimentalists [l| come to such 
a conclusion. Second, one can come to it after the fol- 
lowing consideration. As for excited states produced by 
external pumping, the statement about low density of 
dipoles can be proven by estimates based on Eq. (|A15I) . 
Note also that some part of the excited atoms can be pro- 
duced by collisions in the system. But for a steady-state 
system the density of such atoms can be comparable only 
in the case when there is an equilibrium between the ra- 
diation (photons) and matter (see Ref. [1]). For a BEC 
in dilute vapors of alkali-metals, in which the ultraslow 
light phenomenon is observed, such equilibrium can be 



achieved, e.g., by a special set of mirrors that returns 
photons to the system. In this case, the interaction of 
closely located radiating dipoles must be taken into ac- 
count because it can have a strong impact on the phe- 
nomenon under study. The mentioned effect should be 
taken into account also for rather dense gases in the BEC 
state. 

Of course, all of the estimates introduced in this sec- 
tion are rather approximate. Probably more explicit es- 
timates can be made only from an approach related to 
the Green-function formalism (see, e.g., Refs. @, [I])- 
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